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ABSTRACT 
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Laser  material  damage  experiments  for  this  thesis  were  the  first  ever  conducted  at 
the  new  DoE  Thomas  Jefferson  National  Accelerator  Facility  (TJNAF)  free  electron  laser 
(FEL)  user  laboratory.  In  the  past  only  large-scale  laser  experiments  were  thought  to 
properly  model  weapons  applications.  Experimental  procedures  developed  in  this  thesis 
allowed  a  scaled-down  laser  of  a  few  hundred  Watts  to  characterize  the  damage  from  a 
weapon-scale  one  million  Watt  laser.  The  TJNAF  FEL  has  the  power  of  a  microwave 
oven  concentrated  into  a  beam  the  size  of  a  pencil  lead.  The  unique  TJNAF  FEL  beam 
bombards  the  target  with  a  steady  stream  of  tens  of  millions  of  pulses  per  second  each 
containing  50  million  Watts  of  power  in  a  short  burst  of  4  x  10"13  seconds.  No  other  laser 
combines  these  characteristics,  and  no  experiments  have  previously  been  done  to  explore 
the  effects  of  the  FEL  pulse.  Target  materials  were  obtained  from  the  Naval  Research 
Laboratory  (NRL)  and  from  Naval  Surface  Warfare  Division  (NSWD)  Port  Hueneme. 
Data  were  collected  and  analyzed  using  video  cameras,  optical  microscopes  and  a 
scanning  electron  microscope  (SEM).  This  thesis  has  been  a  productive  cooperation 
among  NPS,  TJNAF,  NRL,  and  NSWD  Port  Hueneme,  to  the  benefit  of  DoD. 
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I.  INTRODUCTION 

This  thesis  argues  that  the  Navy  has  a  need  for  new  weapons  to  provide  ship  self- 
defense  against  sea-skimming  missiles  in  littoral  waters.  Since  the  Navy's  Directed 
Energy  Office  has  already  selected  the  free  electron  laser  (FEL)  for  developmental 
funding,  this  thesis  describes  the  FEL  as  a  candidate  weapon.  The  primary  thrust  of  the 
thesis  is  to  describe  laser  damage  experiments  conducted  at  the  Thomas  Jefferson 
National  Accelerator  Facility  (TJNAF).  These  are  the  first  experimental  tests  that  study 
the  damage  from  a  short  pulsed  laser  at  a  high  repetition  rate  with  a  few  hundred  Watts  of 
average  power.  The  unique  idea  advanced  in  this  thesis  is  that  scaling  rules  can  be 
developed  that  will  allow  these  small-scale  damage  experiments  to  represent  the  damage 
from  a  large,  MW-scale  weapon. 

A.         SHORTCOMINGS  OF  CURRENT  SHIP  SELF-DEFENSE  SYSTEMS 

With  the  collapse  of  the  Soviet  Union  and  the  fading  from  memory  of  the  cold 
war  era  bipolar  world,  a  "new  world  order"  was  supposed  to  have  emerged.  There  were 
rosy  predictions  of  a  world  without  conflict  where  the  nations  of  the  earth  coexisted  in 
peace  and  prosperity.  Now  less  than  ten  years  after  the  fall  of  the  Berlin  Wall,  we  see  that 
those  predictions  were  grossly  optimistic,  and  that  new  challenges  have  arisen  for  the  art 
of  diplomacy  and  for  military  science.  We  have  had  constant  bloodshed  from  the  Gulf 
War  to  the  never-ending  ethnic  violence  in  the  Balkans.  In  addition,  many  countries 
whose  interests  do  not  coincide  with  the  interests  of  the  United  States  seem  to  be 
vigorously  pursuing  programs  to  develop  weapons  of  mass  destruction  and  advanced 


missile  technology.  Many  of  these  countries'  programs  seem  to  have  been  enhanced  by 
the  break-up  of  the  Soviet  Union,  which  has  made  weapons  technology  and  technicians 
available  to  the  highest  bidder. 

The  U.S.  military  by  contrast  has  been  obliged  to  "do  more  with  less"  due  to 
budget  cuts,  while  changing  weapons  procurement  and  tactics  to  meet  emerging  threats. 
For  the  U.S.  Navy,  there  is  no  longer  another  navy  in  the  world  to  pose  a  threat  on  the 
high  seas,  and  the  primary  focus  has  shifted  to  the  littoral  with  emphasis  on  power 
projection  ashore  from  the  sea  and  support  of  land  forces.  This  new  operating 
environment  has  revealed  vulnerabilities  of  U.S.  Naval  forces.  These  vulnerabilities  were 
not  anticipated  when  current  ship's  systems  were  being  designed  and  built  during  the 
cold  war's  uncertain  aftermath.  The  primary  vulnerability  of  U.S.  Navy  ships  is  attack  by 
high-speed  anti-ship  missiles  and  is  exacerbated  by  operating  in  the  littoral  environment. 

1.  Surface-To- Air-Missiles  (SAMs) 

The  first  line-of-defense  for  most  ships  is  the  Standard  Missile,  which  is  launched 
vertically  and  designed  to  destroy  incoming  missiles  at  long  range.  These  systems,  along 
with  ships'  long-range  radars,  are  ideally  suited  to  the  open  ocean  environment  with 
visibility  limited  only  by  the  curvature  of  the  earth.  In  the  littoral  on  the  other  hand, 
conditions  are  often  quite  different.  Because  ships  often  operate  within  a  few  miles  of 
land,  the  detection  range  for  incoming  missiles  is  severely  curtailed  with  subsequent  drop 
in  reaction  time.  Due  to  the  speed  of  modern  anti-ship  missiles,  and  the  short  ranges  of 
the  littoral,  the  Detect-To-Engage  (DTE)  sequence  for  the  Standard  Missile  will  be 
impossible  to  execute  in  the  time  available. 


Furthermore,  these  Standard  Missile  systems  are  not  installed  on  all  Navy  ships, 
but  only  on  the  cruisers  and  the  destroyers,  about  two-thirds  of  the  ships.  The  rest  of  the 
ships  have  no  "long  range"  anti-ship  missile  defense  at  all,  for  either  the  open  ocean  or 
the  littoral.  The  system  is  very  expensive  both  in  terms  of  money  and  in  terms  of  the 
space  required.  Since  a  ship  has  a  very  limited  budget  of  space,  and  the  missile  magazine, 
the  missile  launching  system,  and  associated  radars  are  large  and  heavy,  it  is  impractical 
to  ever  put  them  on  logistics  ships  or  amphibious  ships.  Cruisers  and  destroyers  will  have 
enough  problems  defending  themselves  from  missile  attack  in  the  littoral,  and  will 
therefore  be  completely  incapable  of  defending  other  ships. 

2.         Close-In-Weapon-System  (CIWS) 

All  ships  that  might  operate  in  the  littoral  do  have  the  CIWS,  which  is  a  short- 
range  gun  and  radar  system  designed  to  destroy  an  inbound  missile  just  before  it  gets  to 
the  ship.  CIWS  has  never  been  tested  in  battle  conditions.  Computer  simulations  indicate 
that  it  is  unlikely  the  system  can  destroy  an  anti-ship  missile  at  ranges  greater  than  100- 
200  meters.  Although  the  gun  has  a  nominal  range  of  2000  meters,  penetrator  dispersion 
seriously  degrades  the  effective  range  of  the  system.  CIWS  has  an  average  dispersion  of 
two  to  three  milliradians  or  0.17  degrees,  which  leads  to  a  probability  of  a  bullet  hitting 
the  target  as  shown  in  Figure  (1).  The  system  has  a  probability  of  hitting  the  target  of 
very  nearly  0%  from  its  maximum  range  all  the  way  in  to  about  500  meters.  The 
probability  does  not  rise  above  10%  until  the  inbound  missile  is  within  200  meters  of  the 
ship.  As  the  missile  approaches  the  ship  the  probability  of  a  hit  increases  and  the 
accumulated  hits  can  be  calculated  by  applying  the  appropriate  probabilities  for  each  100- 
meter  increment  of  the  missiles  path. 
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Figure  1.  Single  Round  Intercept  Probability 


Assuming  a  relatively  slow,  non-maneuvering  missile,  and  assuming  it  takes  an 
average  of  eight  hits  to  kill  a  missile,  Figure  (2)  shows  that  the  kill  range  will  be  100-200 
meters  as  mentioned  above.  Killing  the  missile  with  multiple  hits  is  not  the  end  of  the 
story.  Once  the  missile  breaks  up,  the  fragments  can  still  hit  the  ship  and  cause  damage. 
The  computer  simulation  was  extended  to  analyze  the  trajectories  of  post-kill  missile 
fragments. 
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Figure  2.  Cumulative  Hits  vs.  Range 


The  shorter  the  kill  range,  the  more  likely  that  many  missile  fragments  will  strike  the 
ship.  Figure  (3)  shows  this  to  be  the  case,  and  in  fact  at  the  likely  kill  range  for  CIWS, 
about  50%  of  the  missile  fragments  will  hit  the  ship  with  the  average  fragment  having  a 
mass  of  40Kilograms,  a  velocity  of  50  meters/second  and  kinetic  energy  of  50000  Joules. 
If  the  missile  could  be  destroyed  at  1000  meters  or  greater  from  the  ship  the  probability 
of  damage  from  missile  fragments  would  be  small.  CIWS,  however,  is  not  capable  of 
killing  targets  at  those  ranges. 


0.6 


05 


04 


03 


02 


01 


100 


Probability  of  a  Fragment  Striking  the  Ship 
Missile  Speed  600m/s,  Missile  Height  5m,  10000  Fragments  Sampled  Every  100m 


300 


500 


700 


900       1 100 
Range  (m) 


1500 


1700 


Figure  3.  Probability  of  Fragment  Striking  Ship  vs.  Kill  Range 

So  even  under  ideal  conditions  for  CIWS  with  ample  warning  to  commence 
engagement  at  maximum  range,  a  slow,  non-maneuvering  missile  leaves  the  phalanx 
system  inadequate  for  ship  self-defense.  Under  more  realistic  littoral  conditions  with  little 
or  no  advance  warning  of  an  attack,  and  a  fast,  maneuvering  missile  the  CIWS  will  be 
nearly  useless.  A  new  ship  self-defense  system  is  required. 


B. 


DIRECTED  ENERGY  WEAPONS 


One  promising  possibility  is  to  use  beams  of  high  energy  focused  on  an  incoming 
missile  to  destroy  it  at  long  range.  A  major  advantage  of  this  type  of  weapon  is  that  it 
travels  at  300  million  meters  per  second.  Since  the  missile  can  only  travel  at  a  few 


hundred  meters  per  second,  any  maneuver  the  missile  makes  becomes  meaningless. 
There  are  three  major  methods  of  generating  the  large  energies  necessary  to  destroy  an 
anti-ship  missile  in  flight:  Particle  Beams,  High  Powered  Microwaves  (HPM),  and 
Lasers.  All  of  these  technologies  have  been  studied  for  many  decades  to  determine  their 
suitability  as  anti-missile  weapons.  The  advantages  and  disadvantages  of  each  will  be 
discussed  in  turn  to  see  where  the  scientific  research  has  lead  so  far. 

1.  Particle  Beams 

The  idea  of  using  a  high-energy  beam  of  subatomic  particles,  electrons,  protons  or 
neutrons,  to  shoot  down  a  missile  has  been  around  since  the  early  1970's,  but  was  given 
major  consideration  and  funding  as  part  of  the  Strategic  Defense  Initiative  (SDI)  in  the 
1980's.  Of  all  the  possible  Directed  Energy  Weapons  the  Particle  Beam  seems  to  face  the 
most  serious  obstacles. 

The  main  problem  for  an  electron  beam  anti-missile  weapon  is  propagation  loss, 
which  has  three  major  components:  Ionization  Loss,  Bremstrahlung  Loss,  and  Beam 
Divergence  due  to  elastic  scattering.  Ionization  Loss  occurs  when  energy  from  the  beam 
ionizes  surrounding  air  molecules  and  amounts  to  approximately  0.34  MeV/m.  This 
means  that  a  500  MeV  electron  beam  would  travel  less  than  1500  meters  even  if  there 
were  no  other  loss  mechanism.  But  there  are  other  loss  mechanisms.  Bremstrahlung  Loss 
is  caused  when  beam  electrons  are  accelerated  in  curved  paths  around  atomic  nuclei  due 
to  coulomb  interaction.  The  electron  beam  then  loses  energy  by  emitting  Bremstrahlung 
photons.  Bremstrahlung  Loss  alone  would  limit  a  500  MeV  beam  to  a  propagation  path 
of  about  360  meters.  Elastic  scattering  of  the  beam  occurs  when  electrons  "collide"  with 
more  massive  air  molecules  as  in  Bremstrahlung  Loss  and  the  trajectories  of  the  electrons 


are  changed.  This  causes  the  electron  beam  to  spread  as  it  propagates  through  the  air. 
Elastic  scattering  alone  would  cause  a  500  MeV  electron  beam  to  dissipate  in  240  meters 
[Ref.  1].  Actual  experimentation  with  all  loss  mechanisms  working  together  shows  that 
an  electron  beam  only  propagates  a  few  meters  through  the  air. 

2.  High  Powered  Microwaves  (HPM) 

HPM  technology  seems  to  hold  more  promise  than  particle  beams  because  of  the 
maturity  of  the  field  and  we  know  we  can  propagate  them  through  the  air.  But  HPM 
weapons  also  face  challenges.  Assuming  that  an  adequate  source  of  HPM  is  available 
there  are  generally  two  methods  of  employment  to  defeat  incoming  missiles.  One  is  to 
aim  an  intense  beam  at  a  specific  target  to  destroy  it  at  long  range,  and  the  other  is  to 
sweep  a  large  area  in  the  hope  of  disabling  many  targets.  The  second  method  is 
unattractive  for  military  applications  because  of  the  lower  probability  of  kill  for  each 
target  and  because  of  the  increased  likelihood  of  fratricide. 

Once  a  source  is  found  and  a  tactic  is  chosen,  the  HPM  radiation  must  reach  the 
target  by  propagation  through  the  atmosphere.  Atmospheric  attenuation  of  HPM  at  sea  is 
dominated  by  suspended  particles  such  as  dust  and  water  droplets,  which  condense 
during  rain  or  fog.  Figure  (4)  shows  that  in  the  absence  of  fog  or  rain,  a  water  vapor 
absorption  peak  at  22  GHz  and  an  oxygen  absorption  peak  at  60  GHz  dominate 
attenuation  of  HPM.  There  is  a  relative  minimum  in  attenuation  at  35  GHz.  Once  the 
radiation  reaches  the  target  there  are  two  basic  transmission  paths  into  the  electronics  of 
the  missile.  The  "front-door"  path  is  one  designed  to  transmit  microwaves  in  the  normal 
operation  of  the  system,  and  the  "back-door"  path  is  one  of  radiation  coupling  through  a 
path  not  designed  for  transmission.  It  requires  less  energy  on  target  to  defeat  a  missile  by 


the  front-door  path,  but  there  is  certainly  no  guarantee  that  the  enemy  will  cooperate  and 
design  missiles  appropriately. 
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Figure  4.  Atmospheric  Attenuation  of  HPM,  From  Ref.  [2] 


An  additional  problem  for  HPM  weapons  is  the  ease  of  shielding  against  them. 
Herbert,  using  a  digital  watch  and  a  microwave  oven,  conducted  simple  experiments  that 
were  quite  enlightening.  A  three-second  exposure  to  microwaves  caused  an  unshielded 
watch  to  quit  working.  When  a  similar  watch  was  placed  in  a  box  made  of  aluminum  foil 
and  placed  in  the  microwave,  it  was  able  to  withstand  repeated  exposure  including 
individual  exposures  of  as  long  as  90  seconds  with  no  loss  of  function  [Ref  3].  Assuming 
a  relatively  slow  missile  capable  of  300  m/s,  the  shielded  missile  would  travel  27 
Kilometers  during  that  90  seconds  and  the  ship  would  be  dead. 


Another  problem  for  long  wavelength  HPM  is  diffraction.  Beam  spread  angle  is 
proportional  to  wavelength  by  the  relation  6  =  1 .22  X I D,  where  D  is  the  aperture  size. 
Beam  spot  area  is  then  A  =  n  (R0)2  ,  where  R  is  the  range  to  the  target.  With  an  antenna 
diameter  of  10  meters  and  a  wavelength  of  6  cm  (5  GHz),  the  beam  spot  is  4208  m2  at  a 
range  of  5  km.  The  energy  is  spread  over  a  wide  area  and  the  effectiveness  of  the  weapon 
reduced. 

3.  Lasers 

Lasers  are  in  use  in  our  society  today  for  so  many  mundane  purposes  that  we  take 
them  for  granted.  We  use  them  as  pointers  for  presentations  in  offices,  to  play  our  music 
and  movies  at  home,  to  cut  things  in  industry,  and  to  perform  surgery  in  hospitals.  Laser 
weapons  have  excited  the  American  imagination  for  decades.  In  the  1950's  early  science 
fiction  shows  had  "ray  guns",  in  the  1960's  Star  Trek  had  "phasors",  in  the  1970's  and 
1980's  Star  Wars  had  "Death  Stars"  with  lasers  which  could  destroy  entire  planets  at  a 
shot.  The  scientific  community  has  worked  hard  in  conjunction  with  the  military  to 
explore  the  potential  of  the  laser  as  a  weapon.  The  first  proof  of  this  potential  was  the  Air 
Force's  ALL  project  in  the  1970's,  followed  by  the  Navy's  MIRACL  project  in  the  1980's, 
and  the  Air  Force's  ABL  program  in  the  1990's  [Ref.  4].  There  are  many  different  types 
of  lasers  that  have  different  strengths  and  weaknesses,  and  all  must  be  judged  against  the 
stringent  requirements  of  a  shipboard  laser  weapon. 

The  laser  must  be  able  to  emit  enormous  power,  but  be  small  enough  to  fit  on  a 
ship.  It  must  operate  at  a  wavelength  that  propagates  well  through  the  atmosphere,  so  that 
range  does  not  suffer.  It  must  not  produce  dangerous  byproducts  that  cannot  be  disposed 
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of  at  sea.  We  shall  look  more  closely  at  each  one  of  these  criteria  and  see  how  the 
available  lasers  stack  up. 

4.  Atmospheric  Propagation 

It  is  important  to  study  the  transmission  of  laser  radiation  through  the  atmosphere, 
because  certain  wavelengths,  and  therefore  certain  lasers,  can  be  eliminated  on  the  basis 
that  they  cannot  propagate  5  km  to  the  target  with  sufficient  energy  to  destroy  a  target.  A 
non-linear  phenomenon  known  as  thermal  blooming  can  occur  with  high  power  lasers 
when  the  air  is  stagnant  and  the  target  is  approaching  on  a  constant  bearing.  The  air 
between  the  ship  and  the  target  is  heated,  changing  the  index  of  refraction,  and  forming 
an  effective  lens  that  can  disperse  the  laser  beam.  Since  thermal  blooming  is  a  non-linear 
effect,  the  more  power  in  the  beam,  the  faster  thermal  blooming  will  occur.  As  a  result  it 
is  critical  to  find  a  wavelength  with  the  minimum  absorption  and,  therefore,  maximum 
transmission  through  the  atmosphere.  Figure  (5)  shows  atmospheric  transmittance  for 
wavelengths  of  interest  for  lasers,  but  the  finer  detail  is  necessary  to  pick  an  exact 
window  for  optimal  transmission.  There  are  regions  of  high  transmittance  in  the  3-5 
micron  band  and  the  8-12  micron  band,  which  have  been  the  typical  wavelengths  for  IR 
sensors.  There  are  much  smaller  sections  of  the  spectrum,  which  give  close  to  100% 
transmittance  at  1.042  um,  1.06  urn,  1.6  urn,  2.2  urn,  and  3.8um.  Cook  and  Albertine 
have  done  a  more  detailed  study  of  this  issue  in  relation  to  a  maritime  deployed  high- 
energy  laser  weapon  system  (HELWS).  Their  conclusions  are  summarized  in  Figures  (6) 
and  (7),  which  show  the  most  promising  wavelength  to  minimize  absorption  and  the 
thermal  blooming  problem  is  1.042  um. 
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Figure  5.  Atmospheric  Transmittance  for  0-15  |um,  From  Ref.  [5] 
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Figure  6.  Total  Absorption  1/km,  From  Ref.  [6] 
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Figure  7.  Total  Extinction  1/km  From  Ref.  [6] 

The  FEL,  which  can  be  designed  for  a  wide  range  of  wavelengths,  is  the  only 
laser  capable  of  producing  this  exact  wavelength.  In  addition  to  design  flexibility,  once 
the  FEL  is  built  it  can  be  easily  modified  by  changing  the  electron  beam  energy  to 
operate  over  a  range  of  wavelengths  within  a  factor  of  two.  Other  lasers  such  as  chemical 
lasers,  gas  discharge  lasers,  excimer  lasers,  and  x-ray  lasers  are  confined  to  a  specific 
wavelength  or  wavelengths  by  their  generation  mechanism.  So,  an  analysis  of  the  first 
requirement,  atmospheric  propagation,  shows  the  FEL  to  be  a  viable  option  for  a  laser 
weapon  to  defend  against  anti-ship  missiles. 

5.  Size  and  Power 

Experiments  in  conjunction  with  the  MIRACL  program  indicate  that  an  energy 

flux  or  power  density,  O  ~  10  kW/cm2  is  needed  to  destroy  a  missile  with  a  dwell  time  of 
a  few  seconds. 
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°=5=w-  (1) 

where  P  is  power,  A  is  spot  size,  R  is  the  range,  and  6  is  the  beam  half  angle  in  the  far 
field.  So  for  a  spot  size  of  100  cm  ,  a  Megawatt-class  laser  is  required.  By  integrating  O 
with  respect  to  time  we  can  find  the  energy  deposited  at  the  target  in  a  given  time,  or 
fluence,  F. 

F  =  jW  =  Or  ,  (2) 

where  x  is  the  pulse  length  and  <f>  is  assumed  constant  over  the  pulse.  By  substituting  Eq. 
(1)  for  O  and  the  diffraction  limited  0  =  1.22  A./D  into  Eq.  (2)  we  get  an  expression  for 
the  fluence. 

PD\ 

F=021W  (3) 

where  X  is  the  wavelength  and  D  is  the  diameter  of  the  aperture.  Because  of  atmospheric 
extinction  more  energy  is  necessary  at  the  ship  to  obtain  the  required  energy  at  the  target 
[Ref.  7]. 

f=2ixio'j£?t-  (4) 

where  T=  ectR  is  the  atmospheric  transmittance,  and  a  is  the  extinction  coefficient,  F 
[kJ/cm2],  P  [MW],  D  [m],  x  [sec],  Mum],  and  R  [km]. 
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II.  FEL  THEORY 

A.         BASIC  OPERATION 

A  free  electron  laser  (FEL)  takes  advantage  of  a  simple  physical  characteristic. 
When  bunched  electrons  are  accelerated  they  emit  photons  coherently.  A  static  magnetic 
field,  transverse  to  the  direction  of  a  beam  of  relativistic  electrons,  causes  the  electrons  to 
change  directions  or  "wiggle,"  from  side  to  side.  As  the  bunched  electrons  wiggle  they 
give  off  coherent  radiation  with  proper  design  of  the  electron  beam,  the  optical  cavity, 
and  the  magnetic  field.  The  major  components  of  an  FEL  are  a  source  of  relativistic 
electrons,  an  undulator  and  appropriate  optics  for  a  resonant  cavity  [Ref.  8]. 

1.  Electron  Beam 

Electrons  are  produced  by  an  electron  gun,  then  accelerated  to  relativistic 
energies.  The  electron  energy  is,  E  =  ymc  ,  where  y  is  the  Lorentz  factor,  m  is  the 
electron  mass,  and  c  is  the  speed  of  light.  Electron  energies  may  be  from  a  few  MeV  to  a 
few  GeV  depending  on  the  desired  laser  wavelength.  Bunching  of  the  electrons  is 
required  for  optical  gain  to  occur.  The  size  of  the  accelerator  is  a  concern  for  shipboard 
application  of  the  FEL  and  electron  beam  quality  is  also  important.  Electron  beam 
quality  is  a  measure  of  the  spread  in  energies  of  electrons  throughout  the  beam,  and  is 
best  described  (smaller  is  higher  quality)  by  the  ratio  of  the  spread  of  the  Lorentz  factors 
to  the  average  Lorentz  factor,  Ay  /  y.  A  typical  layout  for  an  FEL  with  a  recirculating 
electron  beam  is  shown  in  Figure  (8). 
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Figure  8.  FEL  Electron  Source 


2. 


Undulator 


In  the  undulator,  the  optical  field  is  formed  and  amplified  by  the  interaction  of  the 
relativistic  electron  beam  with  the  spatially  periodic  magnetic  field.  The  magnetic  field  is 
produced  by  the  alternating  polarity  arrangement  of  the  magnets.  A  picture  of  an 
undulator,  sometimes  referred  to  as  a  "wiggler,"  is  provided  in  Figure  (9).  The  undulator 
wavelength,  A-o,  is  the  distance  along  the  beam  axis  between  magnet  pairs,  and  is  given 
by,  A-o  =  L/N,  where  L  is  the  total  length  of  the  undulator  and  N  is  the  number  of 
undulator  periods.  The  magnetic  field  can  be  produced  by  electromagnets,  but  most  FEL 
undulators  are  generally  constructed  with  permanent  magnets  as  in  Figure  (9). 


Figure  9.  FEL  Undulator 
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3.  Optical  Cavity 

The  optical  cavity  is  oriented  along  the  undulator  axis  and  extends  beyond  the 
undulator  length.  In  Figure  (10),  the  arrows  in  the  center  represent  the  undulator  periods, 
which  are  bracketed  by  mirrors  containing  the  optical  field.  During  successive  passes 
through  the  undulator,  the  optical  field  is  amplified  and  some  fraction  of  the  coherent 
radiation  is  allowed  to  escape  and  used  in  the  weapon  system.  The  physics  of  the  FEL 
can  be  described  by  two  equations  -  one  for  the  electron  motion  and  one  for  the  optical 
field  evolution.  These  equations  will  be  described  in  detail  in  the  following  section. 


FREE    ELECTRON  LASER  RESEARCH 

Figure  10.  FEL  Optical  Cavity 

B.         PENDULUM  EQUATION 

It  will  be  shown  that  the  pendulum  equation  describes  the  transfer  of  momentum 
and  energy  between  a  free  electron  and  the  electromagnetic  wave  in  the  magnetic  field  of 
the  undulator.  The  governing  force  equations  for  the  electron  are  the  relativistic  Lorentz 
equations, 

d(yQ)         e     -    — 

-^P  = [£,  +  Px(5r  +  5)],  (5) 

at  mc 

at         mc 

y  "2  =  1  -  p  •  (3  ,  (7) 


17 


where  P  is  the  velocity  of  the  electron  as  a  fraction  of  the  speed  of  light.  The  electrical 
field  and  magnetic  field  of  the  optical  wave  are  given  by, 

Er  =  E0[cosC¥),smC¥),0],  (8) 

Br  =  B0[smm,cosQ¥),0],  (9) 

where  *F  =  kz  -  cot  +  <j> ,  the  wave  number  is  k  =  In/X  ,  co  is  the  angular  frequency,  <j)  is 
the  optical  phase  angle,  and  A,  is  the  optical  wavelength.  The  undulator  static  magnetic 
field  is  represented  by, 

B  =  B0[cos(k0z),sm(k0z),0],  (10) 

where  the  undulator  wave  number  is  k0  =2izl\  ,  and  X^is  the  undulator  wavelength. 
Substituting  the  fields  in  Eq.  (8),  (9),  and  (10)  into  Eq.  (5),  and  resolving  the  transverse 
components  gives, 

"lp  = [£0(l-|3r)(cos(^),-sin(T),0)  +  50P2(-sin(VXcos(VX0)],    (11) 

ot  mc 

where  p\  is  the  transverse  electron  velocity  and  p\,  is  the  electron  velocity  in  the  direction 

of  the  undulator  axis.  For  a  relativistic  electron  E0(\  -  Pr)  «  E0/2y  2 ,  which  can  be 

neglected  for  y  »  1 ,  leaving, 

SfvB  ^  e 

-^2  =  _  — ^p^-sin^oZ^cosC^zXO].  (12) 

ot  mc 

Eq.  (12)  can  be  integrated  in  time  with  the  constant  of  integration  set  to  zero  indicating 
ideal  injection  of  electrons.  The  result  is, 

-         K 

Px  = [cos(£0z),sin(A:0z),0],  (13) 

Y 
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where  we  have  defined  an  undulator  parameter  K  =  eB0\ / 'liunc2 .   The  undulator 
parameter  allows  adjustment  of  the  transverse  electron  velocity  by  adjusting  the 
undulator  magnetic  field  intensity  or  the  spacing  between  magnets. 
Insert  Eq.  (8)  and  Eq.  (13)  into  Eq.  (6)  to  get, 


dy         eE0K 
dt         ymc' 


Y=-7T  = -7^05(^+4)),  (14) 


v  eE  K 

L  =  -V-T<;os(C  +40,  (15) 

y         y  mc 

where  the  electron  phase  is  C,  =  (k  +  k0  )z  -  cat .  The  initial  phase  at  time  t  =  0  is 

C#  =  (k  +  k0  )z0 ,  but  since  k »  k0  we  can  approximate^  «  kz0  =  2%z0 /X ,  which  gives  a 

relation  between  the  electron  phase  and  the  optical  wavelength.  Insert  Eq.  (13)  for  p\ 
into  Eq.  (7)  to  get 

Y"*=l-P-P=l-0x-P,2=Jjfk.  (16) 

Differentiate  both  sides  with  respect  to  time  to  get, 


Y      Y2PsPr 


y  '    l  +  K2  ' 
then  using  Eq.(15)  substitute  P    =  C,  /(k  +  k0)c ,  and  solve  for  C,  , 


(17) 


\  +  K2  (k  +  k0)ceE0K 
P.       y  2fnc 


£,=       2  p°        2°  _cos(^+(|)).  (18) 


Applying  approximations  based  on  k  »  k0 ,  Pz «  1  when  y»  1 ,  and  defining 
dimensionless  time  x  =  ct/L,  we  have  the  pendulum  equation, 

£  =  McosG+4>),  (19) 
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a  = 


4nNeKLE0 

2         2         ' 

y  mc 


(20) 


where  \a\  is  the  dimensionless  optical  field,  N  is  the  number  of  undulator  periods,  L  is  the 

o 

total  length  of  the  undulator,  and  (_)  =  d(J)/dx. 

C.         OPTICAL  WAVE  EQUATION 

The  second  part  of  the  mathematical  description  of  an  FEL  is  the  wave  equation. 
Beginning  with  Maxwell's  wave  equation  for  the  vector  potential  A  such  that 


E=- 


\_dA 
c  dt  ' 


Br  =  VxA, 


V2- 


1   d: 


2    a*  2 


cl  dt 


471    - 

A  =  —J±, 
c 


(21) 
(22) 
(23) 


where  J±  is  the  transverse  current  density  and 


A(z,  t)  =  -E(z,  0[cos(vF)  -  sin(T),0] , 

CO 


(24) 


where  the  optical  phase  is  4*  =  kz  -  cot  +  <j)(z,  t).  Taking  two  spatial  derivatives  and  two 
time  derivatives  of  Eq.  (24)  yields, 


d2A      1  dE 

dz2      k  dz 


k  + 


dz 


[cos(vF),-sin(vF),0] 


ld2ET.  n     1  dE 

—  [sin(4/),cos(xF),0]  + 


k  dz' 


k  dz 


k  + 


d§_ 
dz 


[cos(lF),-sin(xF),0] 


£d2<i>r  n     E 

+  j-^r[cosQ¥)-smC¥)t6\  +  - 


k  + 


d§_ 
dz 


[-sin(^),-cos(T),0],  (25) 
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and 


co  d2A      dE 
c  dt2  "  dt 


dt 


-co 


d2E 


[cos(vF),-sin(vF),0]  +  — —  [sinOF),cos(¥),0] 


dt' 


+ 


dE 
dt 


cfy 

dt 


-co 


sf 


[cos(¥)  ,-sin(¥),0]  +  E-^y[cosC¥) ,-sin(¥),0] 


+  E 


3<J> 

dt 


-co 


[-sin(vP),-cos(T),0]. 


(26) 


Next,  assume  the  optical  phases  and  amplitudes  are  slowly  varying  in  time  and 
space  such  that  dE/dz  «  kE ,  dty/dz  «  k§  ,  dE/dt  «  (oE ,  dty/dt  «  cocb ,  and 
co  =  kc.  Eliminate  second-order  terms  and  insert  Eq.  (25)  and  Eq.  (26)  into  Eq.  (23)  to 
get 


V2- 


1    d- 


2    a*2 


cl  dt 


A  *  2 


+  2E 


dE     \dE_ 
dz      c  dt 


[cos(vP),sin(vF),0] 


cfy      ]_d§_ 
dz     c  dt 


Arrr 

f-sin(vF)-cos(4/),0l  *  -  —  J,. 
J  c 


(27) 


The  current  density  for  a  single  electron  is  J±  =  -ec$L.  Substituting  Eq.  (13)  for  (3X  and 
introducing  new  time  and  space  coordinates,  x  =  ct/L,  and  z  =  z-ct,  we  have, 


\_dE_ 
Ldx 


[cos(vF),sin(4'),0]  +  2£  7^7  [-sin(vF),-cos(lP),0] 


4neK 

Y 


[cos(£0z),sin(£0z)J.  (28) 


The  dimensionless  time  t  progresses  from  zero  to  one  from  the  beginning  of  the 
undulator  to  the  end,  and  the  coordinate  z  follows  a  point  on  the  optical  pulse  traveling 
at  speed  c.  Insert  Eq.  (28)  into  Eq.  (27),  take  the  average  (...)  over  many  sample 
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electrons,  and  multiply  by  the  electron  particle  density,  p,  giving  the  equations  of  motion 
for  the  optical  field  amplitude  and  phase, 

dE        InepLK 


dx  y 

d<f)      2nepLK 


cos 


£+♦)),  (29) 


sin 


£+♦)),  (30) 


dx  y 

where  C,  =  (k  +  k0)z  -  cor  is  the  electron  phase.  Use  the  complex  electric  field  E  =  E0e®  , 
to  get 

dE        InepLK  r,  \      ,  .  vl  ,.         InepLK  ,  iC 


[{cosK+d,)>-/(sine+(t»)»'*  =-2=2±{e<).  (31) 


dx  y  y 

Substituting  the  definition  of  the  dimensionless  optical  field  a  =  \a\e*  from  Eq.  (13)  we 
get  the  final  form  of  the  wave  equation 

da      °  /    ,r\ 

_  =  a  =  _y(^)5  (32) 

where  the  dimensionless  current  is  j  =  8n  2e2pNK2L2/y  3mc2.  The  growth  of  the  optical 
field,  and  therefore  the  development  of  laser  energy,  is  dependent  on  the  dimensionless 
current  and  the  average  electron  phase.  So  we  need  both  a  current  j  and  electron 

bunching  (e~'? )  *  0  to  produce  laser  energy. 
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III.  FEL  DESIGN  AND  PARAMETERS 

A.         TJNAF  FEL 

The  most  powerful  FEL  ever  operated  is  at  the  Thomas  Jefferson  National 
Accelerator  Facility  (TJNAF).  In  April  1999,  the  laser  operated  continuously  at  710 
Watts  of  average  power  for  six  hours.  The  FEL  is  designed  to  achieve  1  kW  and  will 
probably  soon  do  so.  Near  term  modifications  now  in  the  planning  stage  will  boost  the 
power  to  20  kW.  Table  (1)  shows  the  parameters  of  the  TJNAF  FEL  and  compares  them 
to  the  requirements  for  a  shipboard  anti-missile  defense  weapon  [Ref.  9]. 


Parameter 

TJNAF  FEL 

Weapon  FEL 

Average  Power 

P=1KW 

P=1MW 

Average  Current 

7=5mA 

1  =  900  mA 

Electron  Energy 

ymc2  =  48  MeV 

ymc2  =  100  MeV 

Lorentz  Factor 

y  =94 

7  =  196 

Undulator  Field 

£0=.38T 

50=.27T 

Undulator  Wavelength 

\)  =  2.7  cm 

\,  =  4cm 

Undulator  Periods 

N=40 

N=25 

Undulator  Length 

L=  108  cm 

I  =100  cm 

Undulator  Parameter 

£=0.96 

K=\ 

Rayleigh  Length 

z0  =  40  cm 

z0=  10  cm 

Electron  Charge/Bunch 

7e/c=60pC 

7e/c=1800pC 

Peak  Current 

7  =  60  A 

7=  600  A 

Electron  Beam  Radius 

rb=  1 00  urn 

rb=  300  urn 

Pulse  Length 

x  =  0.4  ps 

x  =  3  ps 

Pulse  Repetition  Rate 

PRR  =  18.7/37.4/74.85  MHz 

PRR  =  500  MHz 

Optical  Resonator  Length 

Lopt  =  8.01m 

Lopt=  10  m 

Output  Coupling 

10% 

10% 

Resonator  Cavity  Losses 

Q  =  <0.5%  /  pass 

Q=  <0.5%/pass 

Optical  Wavelength 

X  =  3-6  jam 

X=  1  urn 

Table  1.  Comparison  of  TJNAF  FEL  with  a  weapon  FEL 
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The  big  differences  are  increases  in  the  (i)  peak  current  by  a  factor  of  1 0,  (ii)  the 
repetition  rate  by  a  factor  of  7,  (iii)  the  electron  beam  energy  by  a  factor  of  2,  and  (iv)  the 
pulse  length  by  a  factor  of  7. 

B.         MODELING  A  MW  FEL 

As  we  can  see  from  Table  (1)  there  is  still  a  long  way  to  go  from  the  current  FEL 
to  a  weapons  application  FEL.  Since  we  do  not  have  a  MW-class  FEL  to  perform 
experiments,  we  must  turn  to  scaling  to  determine  the  effectiveness  of  a  FEL  weapon 
without  spending  hundreds  of  millions  of  dollars.  As  mentioned  in  Chapter  I,  power 
density  <E> «  10  kW/cm  over  a  spot  size  A  «  100  cm  is  required  to  destroy  a  missile  with 
a  dwell  time  of  a  few  seconds.  By  scaling  we  can  get  an  appropriate  spot  size  for  any 
laser  power.  With  a  power  density  of  O  =  10  kW/cm  ,  a  100  W  laser  must  use  a  spot  size 
of  1  mm  ,  and  a  1  kW  laser  must  use  a  spot  size  of  1 0  mm  . 

Scaling  will  only  work,  however,  if  the  heat  diffusion  is  independent  of  spot  size. 
Schriempf  gives  a  detailed  treatment  of  diffusion,  but  a  summary  is  useful.  The  heat  flow 
equation  is 

pC—  =  KV2T  +  A  (33) 

where  p  is  the  material  density,  C  is  the  specific  heat,  T  temperature,  t  time,  K  thermal 
conductivity  and  A  is  the  heat  produced  per  unit  volume  per  unit  time  [Ref.  10].  For  a 
semi-infinite  solid  with  no  phase  change  and  material  properties  that  are  temperature- 
independent,  Eq.  (33)  reduces  to 

2^-121-^  (34) 

dz2     k  dt        K' 
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where  the  thermal  diffusivity  is  k  =  K/pC.  The  heat  produced  per  unit  volume  per  unit 
time  is  given  by 

A(z,t)  =  (1  -  tyQityxe-"2 ,  (35) 

where  9?  is  the  reflectivity,  O(t)  is  the  power  density  at  depth  z,  and  a  is  the  absorption 
coefficient.  At  optical  wavelengths,  for  most  solid  materials,  the  absorption  coefficient  is 
very  large  a  «  1 06/cm.  The  exponential  term  approaches  zero,  and  A  is  small  compared 
to  the  time  and  space  derivative  terms.  As  a  result  we  can  assume  A  «  0  and  Eq.  (35) 
becomes 

d2T      1  dT 


—  —  =  0.  (36) 


The  solution  to  Eq.  (36)  is 


dzl     k  dt 


T(z,t)  =  -^ierf 


z 


(37) 


where  Oo  is  the  constant  power  density  on  the  target  and  the  characteristic  thermal 

diffusion  length  is  D  =  2vk7  which  represents  the  distance  required  for  7  to  drop  to  Me 
times  its  initial  value.  In  the  semi-infinite  approximation  we  ignored  radial  heat  flow,  but 
in  order  for  this  to  be  valid  the  spot  size  must  be  much  larger  than  D,  or  the  target 
diameter  d  <  D.  If  these  conditions  are  not  met,  we  expect  heat  to  diffuse  outside  of  the 
laser  spot  and  the  spot  will  not  be  heated  effectively. 

One  of  the  targets  to  be  irradiated  was  Al-6061 .  A  calculation  of  the  thermal 
diffusion  length  associated  with  heating  the  sample  to  melting  temperature  was 
performed  with  p  =  2700  Kg/m3,  C  =  896  J/Kg-K,  K  =  180  W/m-K,  k  =  7.44  x  10"5  m2/s, 
Tm=  855K,  Oo  =  108  W/m2  and 
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D  =  2Jkx,  (38) 


7iK2Ar2 

t=    ,_2      ,  (39) 


Where  t  is  the  time  required  to  bring  the  material  from  ambient  temperature  to  melting 
temperature  and  AT"  is  the  temperature  change.  The  result  for  Al-6061  was  D  =  1.79  mm. 
In  order  to  melt  through  an  aluminum  sample,  the  scaled  laser  spot  must  have  an  area 
A  »  ti(1  .79  mm)  =10  mm  or  the  target  itself  made  small  with  d  «  D .  These 
calculations  were  experimentally  verified  with  samples  of  Al-6061.  Using  d  =  1  cm,  the 
sample  was  irradiated  with  no  melting  after  several  minutes.  Other  samples  with 
d  =  2  mm  were  melted  in  a  few  seconds. 

Another  target  to  be  irradiated  was  Slip-cast  Fused  Silica  (SiC^).  A  calculation  of 
the  thermal  diffusion  length  associated  with  heating  the  sample  to  its  melting  temperature 
was  performed  with  p  =  2200  Kg/m3,  C  =  920  J/Kg-K,  K  =  1 .26  W/m-K,  k  =  5  x  1 0"7 
m2/s,  Tm=  1980K  [Ref.  11],  O0  =  108  W/m2  and  the  result  was  D  =  0.021  mm.  Therefore, 
with  the  insulating  material,  fused  silica,  the  scaled  laser  spot  must  have  an  area 
A  »  7i(0.021  mm)2  =  0.0014  mm2,  which  is  attained  with  the  1  mm2  beam  used  at 
TJNAF. 

C.         PULSE  TRAIN 

Since  there  is  growing  interest  in  lasers  for  many  scientific  and  engineering 
applications,  research  into  short  pulse  effects  has  increased  sharply  in  the  last  few  years. 
The  pulse  train  of  an  FEL  is  different  from  any  other  laser  and  its  interaction  with  matter 
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at  high-average  power  has  not  been  studied.  The  TJNAF  FEL  has  a  pulse  length  of 
t  =  0.4  ps  and  a  repetition  period  of  T  =  27  ns  illustrated  in  Figure  (11).  The  duty  cycle 
D  is  the  fraction  of  time  the  laser  is  actually  irradiating  the  target, 


x       4xl0"13s  , 

D  =  -  =  ——-r-  =  L5x\o-5. 

T     2.7  x  10    s 


The  peak  power  in  each  micropulse  P  is 


«      P        HOW 
P  =  -=  .=47MT, 

D     1.5  xlO-5 


Where  P  is  the  average  power  of  710  W. 


(40) 


(41) 


m 


m 


Figure  11.  FEL  Pulse  Format 

Comparing  the  TJNAF  FEL  to  another  short  pulse  laser  is  instructive.  The 
Lawrence  Livermore  National  Labs  (LLNL)  1.053-um  Ti:sapphire  CPA  system  [Ref.  12] 
has  a  pulse  length  x  =  0.4  ps,  but  a  pulse  repetition  rate  of  only  10  Hz,  so  the  period  is 
T  =  0. 1  s,  and  the  peak  power  is  P  =  2.5TW  .  The  duty  factor  is 


^     4x10  13  s     A      ■    12 

D  = =  4xl0"12 

01s 


(42) 


so  that  the  average  power  is 
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P  =  (2.5  x  1012^)(4  x  1(T12)  =  \0W .  (43) 

Note  that  the  LLNL  laser  has  a  much  higher  peak  power  than  the  TJNAF  FEL,  but  the 
TJNAF  FEL  has  more  than  seventy  times  the  average  power  because  of  its  high  duty 
cycle.  The  experiments  detailed  in  the  following  section  were  conducted  to  study  the 
effects  of  the  unique  FEL  pulse  format  in  laser-matter  interaction  on  small  samples. 
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IV.  FEL  EXPERIMENTS 

A.         LABORATORY  LAYOUT 

Experiments  for  this  thesis  were  conducted  at  the  Thomas  Jefferson  National 
Accelerator  Facility  (TJNAF),  which  is  operated  by  the  Southeastern  Universities 
Research  Association  (SURA)  for  the  U.  S.  Department  of  Energy  (DOE).  Once  a  laser 
beam  has  been  created,  as  described  in  chapter  two,  the  beam  is  sent  by  a  low-loss  optical 
path  to  user  laboratories  in  another  part  of  the  building.  User  laboratory  number  one  was 
used  for  all  experiments  described  below.  M.  Shinn,  S.  Benson,  B.  Yunn,  G.  Neil,  K. 
Jordan  and  J.  Gubeli  operated  the  laser  and  other  laboratory  equipment.  The  optical 
bench  was  configured  with  a  focusing  lens,  a  sample  holder  and  two  video  cameras.  One 
camera  was  focused  on  the  front  of  the  sample  and  one  on  the  back  of  the  sample  as 
shown  in  Figures  (12)  and  (13). 
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Figure  12.  Back  view  of  sample  in  user  lab  arrangement 
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B. 


Figure  13  Front  view  of  sample  in  user  lab  arrangement 


EXPERIMENTAL  PROCEDURE 


Samples  were  irradiated  by  a  laser  beam  with  wavelength  X  =  4.825  urn  through  a 
calcium  fluoride  lens  with  a  focal  length  of  300  mm.  The  pulse  repetition  frequency 
(PRF)  was  37.425  MHz  and  the  average  power  meter  in  the  user  lab  read  100-103  W 
with  an  error  of +/-  5  W.  Since  a  lens  focused  the  beam,  the  beam  area  decreased  with 
distance  along  the  direction  of  propagation  to  a  minimum  waist  radius  of  wo  =  80  jam  at 
the  focal  point.  Paraxia,  a  beam  propagation  code,  was  used  to  model  the  beam 
diffraction  and  find  the  target  position  giving  the  desired  intensity  of  10  kW/cm  [Ref. 
13].  Figure  (14)  shows  a  graph  of  irradiance  versus  distance  from  the  focal  point,  with 
the  negative  numbers  indicating  positions  in  front  of  the  focus. 
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Figure  14.  Irradiance  vs.  Distance  from  Focal  Point,  From  Ref.  [13] 

An  irradiance  of  10  kW/cm  occurs  when  the  sample  position  is  between  25-30  mm  in 
front  of  the  focal  point.  The  samples  were  actually  placed  26  mm  in  front  of  the  focus. 
The  burn-through  time  was  determined  by  observing  a  signal  on  a  power  meter  placed 
behind  the  samples,  and  by  watching  for  the  presence  of  coherent  harmonics  in  the 
visible  spectrum  on  an  iris  placed  approximately  15  cm  behind  the  samples. 


C. 


DESCRIPTION  OF  RESULTS 


The  Naval  Research  Laboratory  (NRL)  provided  the  sample  materials,  which 
included  Slip-cast  Fused  Silica  (SiC^),  Polyimide  Fiberglass,  and  F2  Epoxy.  The  last  two 
were  from  the  same  batch  as  the  material  used  for  previous  continuous  wave  (C  W)  laser 
penetration  studies. 
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1.  Sample  #1  Slip-cast  Fused  Silica 

The  sample  provided  by  NRL  is  6.9  cm  by  7.4  cm  and  varies  in  thickness  from 
0.9  cm  to  1.9  cm.  The  front  of  the  sample  after  irradiation  is  shown  in  Figure  (15). 


Figure  15.  Slip-cast  fused  silica  sample 

The  sample  was  irradiated  eight  times  and  the  results  are  summarized  in  Table  (2). 


Run 
Number 

Average 
Intensity  (kW/cm2) 

Exposure 
Time  (s) 

Penetration  Depth 
Rate  (mm/s) 

1 

9 

9 

0.20 

2 

9 

110 

0.081 

3 

9 

13 

0.20 

4 

9 

24 

0.125 

5 

9 

41 

0.081 

6 

9 

2 

0.35 

7 

500 

3 

3.0 

8 

500 

11 

3.0 

Table  2.  Irradiation  of  Slip-cast  Fused  Silica,  After  Ref.  [14] 
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The  irradiations  were  done  from  left  to  right.  The  last  two  runs  were  done  with 
the  sample  at  the  focus  instead  of  26  mm  in  front  so  that  the  beam  waste  radius  was  w  = 
wo  =  80  urn  and  the  intensity  was  500  kW/cm  .  The  last  two  runs  were  conducted  to 
investigate  the  effects  of  much  higher  power  density.  As  shown  in  Figure  (15),  the  first 
six  runs  were  along  the  top  of  the  sample  and  the  last  two  were  approximately  1  cm 
below.  The  second,  seventh  and  eighth  runs  penetrated  the  entire  0.9-cm  thickness  of  the 
Fused  Silica  material.  Figure  (16)  shows  the  effects  of  exposure  time  on  penetration 
depth  rate. 
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Figure  16.  Exposure  Time  vs.  Penetration  Depth  Rate  for  Fused  Silica 
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As  the  exposure  time  was  increased  more  smoke  and  debris  filled  the  hole  blocking  the 
path  of  the  laser  beam  and  causing  the  penetration  depth  rate  to  decline  over  time.  In  the 
future,  we  can  explore  whether  altering  the  FEL  wavelength  during  sample  irradiation 
improves  penetration  depth  rates  through  smoke  and  debris. 

Since  the  second  run  was  the  only  one  to  punch  completely  through  the  material 
that  was  also  conducted  at  the  primary  power  density  of  interest,  we  shall  take  a  closer 
look  at  it.  A  digital  picture  of  run  two  was  taken  through  an  optical  microscope  as  shown 
in  Figure  (17). 


1  mm 


Figure  17.  Close-up  of  damage  to  slip-cast  fused  silica  in  run  2 

Although  the  beam  diameter  was  only  1 . 1  mm,  the  melted  portion  at  the  surface  of  the 
sample  measured  5  mm  in  diameter.  The  hole  is  tapered  with  the  melted  portion  on  the 
back  of  the  sample  measuring  only  2  mm  in  diameter.  The  reasons  for  this  are  twofold. 
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First,  the  beam  profile  is  Gaussian  in  nature  with  the  highest  intensity  in  the  center  of  the 
beam  and  intensity  down  by  1/e  at  the  beam  radius.  Second,  the  sample  face  is  26  mm  in 
front  of  the  focus  and  the  sample  back  is  17mm  from  the  focus,  so  the  beam  size  is 
decreasing  as  it  proceeds  through  the  material  as  shown  in  Figure  (18). 
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Figure  18.  Sketch  of  beam  focusing  effect 

The  vertical  scale  in  Figure  (18)  is  exaggerated  by  a  factor  of  five  with  respect  to  the 
horizontal  scale  in  order  to  demonstrate  the  effect. 

Examination  of  the  hole  from  run  two  through  an  optical  microscope  reveals  a  1- 
mm  thick  layer  of  melted,  and  rehardened,  SiC>2  filling  the  hole  at  the  back  of  the  sample. 
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It  was  clear  from  the  video  and  the  rear  power  meter  that  burn-through  occurred  in  run 
two,  but  melted  material  solidified  and  sealed  the  hole  at  the  back  of  the  sample.  A 
picture  of  the  back  of  the  target  taken  through  a  Scanning  Electron  Microscope  (SEM), 
Figure  (19),  shows  the  hole  from  run  seven  is  fairly  irregular  with  a  great  deal  of  debris. 


Figure  19.  SEM  photograph  of  damage  to  Fused  Silica  sample  from  run  7 

The  volume  of  the  hole  in  run  two  is  estimated  by, 

V=  ^mmn[R(z)]2dz,  (44) 

where  the  radius  approximately  changes  linearly  as 

R{z)  =  0.53mm  -0.02z,  (45) 

which  gives  a  volume  of  V-  5.6  mm3.  By  a  similar  calculation,  the  volume  of  the  entire 
damaged  region,  including  the  melted  and  rehardened  portion,  is  estimated  to  be 
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V—  92  mm  .  Based  on  the  density  of  fused  silica  of  p  =  2.2  g/cm  ,  the  amount  of 
material  removed  was  0.012  g,  and  the  amount  of  material  damaged  was  0.20  g.  The  heat 
energy  deposited  during  run  two  is  given  by, 

E  =  Pt  =  <S>Ax  ,  (46) 

which  gives  E  =  9.7  kJ  deposited  during  the  1 10  second  run.  The  heat  of  ablation  is  then 
48  kJ/g  for  the  damaged  area.  In  future  studies  this  value  should  be  correlated  to  the 
values  obtained  for  similar  materials  in  much  larger  scale  CW  laser  experiments  in  the 
past,  with  appropriate  compensation  for  experimental  environment  and  wavelength. 

2.         Sample  #2  Polyimide  Fiberglass 

The  sample  provided  by  NRL  is  1 1.4  cm  by  10.1  cm  and  2  mm  thick.  The 
damaged  area  of  the  sample,  after  irradiation,  is  shown  in  Figure  (20). 


Figure  20.  Polyimide  Fiberglass  target 

The  sample  was  irradiated  three  times  and  the  results  are  summarized  in  Table  (3). 
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Run 

Number 

Average 
Intensity  (kW/cm2) 

Exposure 
Time  (s) 

Penetration  Depth 
Rate  (ram/s) 

1 

9 

7 

0.28 

2 

9 

2 

0.90 

3 

9 

1 

1.1 

Table  3.  Irradiation  of  Polyimide  Fiberglass,  After  Ref.  [14] 

The  irradiations  were  done  from  left  to  right  with  the  sample  26  mm  in  front  of 
the  focus  of  the  beam.  Only  the  first  run  achieved  burn-through  of  the  material,  with  the 
entry  hole  3  mm  in  diameter  and  the  exit  hole  1 .5  mm.  All  three  holes  show  significant 
charring  adds  an  additional  term  to  the  heat  transport  equation  and  impedes  ablation. 
Investigation  with  an  optical  microscope  reveals  a  raised  lip  of  material  around  the  face 
of  the  hole  that  does  not  appear  on  the  fused  silica  sample  and  much  more  roughness  as 
observed  in  Figure  (21)-(23).  The  charred  region  extends  to  a  diameter  of  8.0  mm  for  run 
one,  6.5  mm  for  run  two,  and  5.4  mm  for  run  three.  The  lip  height  is  0.3  mm  for  run  one, 
0.1  mm  for  run  two,  and  0.05  mm  for  run  three.  These  measurements  indicate  that  as  the 
dwell  time  increases,  the  radial  extent  of  the  damage  area  increases,  and  more  material  is 
deposited  around  the  edge  of  the  hole.  There  is  no  evidence  of  melted  and  rehardened 
material  present  in  the  holes  as  found  with  the  fused  silica  indicating  a  different 
mechanism  for  ablation  in  the  two  samples. 
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Figure  21.  Close-up  of  damage  to  Polyimide  fiberglass  in  Run  1 


Figure  22.  Close-up  of  damage  to  Polyimide  fiberglass  in  Run  2 
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Figure  23.  Close-up  of  damage  to  Polyimide  Fiberglass  in  Run  3 

3.  Sample  #3  F2  Epoxy 

The  sample  provided  by  NRL  is  10.0  cm  by  1 1 .5  cm  and  1 .5-mm  thick  including 
a  1.6-cm  thick  polyurethane  foam  backing.  The  damaged  area  of  the  sample,  after 
irradiation,  is  shown  in  Figure  (22). 


Figure  24.  F2  Epoxy  sample 

The  sample  was  irradiated  three  times  and  the  results  are  summarized  in  Table  (4). 
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Run 

Number 

Average 
Intensity  (kW/cm2) 

Exposure 
Time  (s) 

Penetration  Depth 
Rate  (mm/s) 

1 

9 

11 

0.10 

2 

9 

6 

0.12 

3 

9 

3 

0.10 

Table  4.  Irradiation  of  F2  Epoxy,  After  Ref.  [14] 

In  each  case,  it  appears  that  the  F2  Epoxy  was  completely  penetrated  and  the 
ablation  of  the  foam  backing  had  begun,  but  not  completed.  The  videotape  showed 
flames  engulfing  the  upper  portion  of  the  sample  and  Figure  (24)  shows  the  black  charred 
area  extending  to  the  edge  of  the  sample.  Significant  charring  was  evident  when  the 
sample  was  viewed  with  the  optical  microscope,  very  similar  to  the  Polyimide  sample. 
There  was  also  evidence  of  some  melting,  but  not  as  much  as  occurred  in  the  Fused  Silica 
sample.  The  holes  appear  to  be  filled  with  the  charred  debris  of  the  polyurethane 
backing,  making  hole  depth  measurements  difficult  and  rendering  penetration  depth  rates 
unreliable.  Figures  (25)-(27)  show  the  details  of  runs  one  through  three. 


Figure  25.  Close-up  of  damage  to  F2  Epoxy  in  run  1 
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Figure  26.  Close-up  of  damage  to  F2  Epoxy  in  run  2 


Figure  27.  Close-up  of  damage  to  F2  Epoxy  in  run  3 

The  damaged  region  extends  to  a  diameter  of  1 1 .3  mm  for  run  one,  7.5  mm  for  run  two, 
and  5.2  mm  for  run  three.  There  is  a  lip  around  each  of  the  holes,  but  much  smaller  than 
the  polyimide  sample  showed.  The  lip  for  run  one  was  0.05  mm  and  for  runs  two  and 
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three  the  lip  was  too  small  to  measure  with  the  optical  microscope  mechanism.  These 
measurements  indicate  that  as  dwell  time  increases  the  radial  extent  of  the  damage  area 
increases,  and  more  material  is  deposited  around  the  edge  of  the  hole.  When  wind  is 
added  to  the  test,  we  may  observe  that  debris  in  the  hole  is  removed  during  irradiation. 
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V.  CONCLUSIONS 

In  this  thesis,  it  was  established  that  a  new  self-defense  weapon  is  needed  to 
counter  current  and  future  threats  to  U.  S.  Navy  ships.  The  FEL  was  proposed  as  a 
possible  alternative  to  the  inadequate  CI  WS  because  of  its  speed-of-light  "bullets", 
unlimited  magazine,  tunability,  and  power.  The  basic  physics  of  the  FEL  was  described 
through  the  derivations  of  the  pendulum  equation  and  the  optical  wave  equation.  The 
Navy's  Directed  Energy  Office  has  already  selected  the  FEL  for  developmental  funding, 
so  experimentation  is  required  to  determine  the  direction  development  should  take.  This 
thesis  describes  the  first  measurements  of  laser  damage  from  the  newly  developed 
TJNAF  FEL  and  the  results  could  provide  the  basis  for  new  directions  for  directed  energy 
weapons  design. 

A.         SCALING 

The  TJNAF  FEL,  which  is  capable  of  several  hundred  Watts  of  continuous 
average  power,  was  used  to  simulate  the  damage  from  a  MW-class  weapon  by  focusing 
the  beam  to  a  smaller  spot  size.  The  eventual  goal  is  to  develop  scaling  rules  that  will 
reliably  predict  the  damage  of  a  large  laser  without  having  to  bare  the  enormous  cost  of 
building  the  large  laser  first.  The  experimental  data  shows  that  the  scaling  concept  with 
thermal  diffusion  calculations  is  promising.  More  detailed  experiments  varying 
wavelength,  power,  and  spot  size  may  be  able  to  produce  scaling  laws,  which  would  be 
invaluable,  for  future  weapons  designers. 
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B.  FEL  PULSE  FORMAT 

The  extremely  short  sub-picosecond  pulse  length  of  the  FEL  beam  is  a  result  of 
the  electron  bunches  described  in  the  section  on  FEL  background  physics.  The  TJNAF 
FEL  has  a  unique  pulse  format  with  a  rapid  sequence  of  short,  powerful  pulses.  The  peak 
power  in  each  pulse  is  about  50  MW  lasting  for  only  about  one-half  picosecond  coming 
at  a  rate  of  37  MHz.  Other  studies  have  shown  that  such  short  pulses  may  give  as  much 
as  a  factor  often  advantage  in  reduced  fluence  required  to  produce  damage  [Ref.  15]. 
The  experiments  conducted  for  this  thesis  began  to  collect  data  to  show  whether  this 
advantage  exists,  but  further  experimentation  will  be  required. 

C.  FUTURE  EXPERIMENTS 

The  TJNAF  FEL  is  scheduled  for  an  upgrade  to  20  kW  of  power  that  will  allow 
more  flexibility  in  scaling  experiments  and  further  tests  of  scaling  itself.  Additional  plans 
are  for  experiments,  which  include  wind  passing  over  the  samples,  weighing  of  the 
samples  before,  and  after  each  run,  new  wavelengths,  changing  wavelength  during 
irradiation,  new  pulse  formats,  and  other  sample  materials.  As  experimental  procedures 
are  refined  and  the  amount  of  data  increases,  more  thorough  analysis  of  the  FEL  beam 
and  comparison  to  other  lasers  will  become  possible. 
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